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ABSTRACT 

I 

Employing  a  technique  that  eliminates  external  field  sources  and  the 
effects  of  strike  aliasing,  we  have  extracted  from  marine  survey  data  the 
intermediate  wavelength  magnetic  anomaly  field  for  H  in  the  North 
Pacific.  A  strong  correlation  exists  between  this  field  and  the  MAGSAT  field 
although  we  can  detect  a  directional  sensitivity  in  the  MAGSAT  field.  The 
intermediate  wavelength  field  is  correlated  to  tectonic  features.  Island  arcs 
appear  as  positive  anomalies  of  induced  origin  likely  due  to  variations  in 
crustal  thickness.  Seamount  chains  and  oceanic  plateaus  also  are  manifested 
by  strong  anomalies.  The  primary  contribution  to  many  of  these  anomalies 
appears  to  be  due  to  a  remanent  magnetization. 

The  source  parameters  for  the  remainder  of  these  features  are  presently 
unidentified  ambiguous.  This  study  indicates  that  the  sea  surface  field  is  a 
valuable  source  of  information  for  secular  variation  analysis  and  the 
resolution  of  intermediate  wavelength  source  parameters. 
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THE  INTERMEDIATE  WAVELENGTH  MAGNETIC  ANOMALY  FIELD 
OF  THE  NORTH  PACIFIC  AND  POSSIBLE  SOURCE  DISTRIBUTIONS 

Introduction 

The  object  of  this  study  is  to  investigate  the  feasibility  of  extracting 
intermediate  wavelength  magnetic  anomalies  of  wavelengths  between  3000  km  and 
300  km  from  total  field  data  acquired  in  marine  magnetic  surveys.  The  long 
term  objectives  of  the  research  are  to: 

(1)  determine  the  distribution  of  intermediate  wavelength  anomalies  over 
the  world's  oceans, 

(2)  determine  the  extent  to  which  MAGSAT  describes  the  distribution  of 
these  anomalies, 

(3)  determine  the  sources  of  the  intermediate  wavelength  magnetic  field. 

Examination  of  the  intermediate  wavelength  magnetic  anomaly  field  of  the 

ocean  basins  offers  the  possibility  of  constraining  some  of  the  ambiguities 
which  affect  magnetic  anomaly  analysis  over  continental  areas.  For  the  most 
part  present  models  for  the  development  of  oceanic  lithosphere  are  relatively 
simple.  The  age  of  the  oceanic  lithosphere  is  well  known  and  models  for  the 
lithospheric  structure  are  simple  and  supported  by  an  extensive  data  set. 
Thermal  sources  and  gradients  have  been  measured  and  delineated.  Furthermore 
isostatic  constraints  require  a  bathymetric  expression  for  bulk  changes  in 
lithospheric  structure.  Therefore,  the  variables  which  might  determine  the 
magnetic  source  distribution  are  relatively  well  constrainted  in  oceanic 
lithosphere. 

In  the  following  pages  we  will  describe  the  techniques  which  were  devel¬ 
oped  to  recover  the  intermediate  wavelength  total  field  anomalies  over  the 
North  Pacific  from  marine  survey  data.  We  feel  that  the  investigation  has 
been  extremely  successful  though  much  remains  to  be  done  in  extending  the 
areas  of  the  study  to  other  oceanic  basins  and  tectonic  features.  Also,  the 
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sources  of  the  observed  intermediate  wavelength  anomalies  mist  be  studied  in 
far  greater  detail  than  has  been  possible  with  the  limited  resources  available 
in  this  pilot  study.  The  technique  is  capable  of  successfully  separating 
intermediate  wavelength  total  field  anomalies  and  of  providing  a  high  resolu¬ 
tion  record  for  secular  variation  during  the  period  1960  to  1980.  Further¬ 
more,  the  sea  surface  data  set  provides  a  higher  resolution  than  the  MAGSAT 
field,  due  to  the  closer  proximity  of  the  sea  surface  data  to  the  lithospheric 
source  bodies  and  eliminates  the  directional  or  longitudinal  attenuation  which 
is  characteristic  of  present  processing  for  the  polar  orbiting  MAGSAT  data. 
The  results  of  the  North  Pacific  study  show  that  the  anomalies  observed  in  the 
MAGSAT  and  P0G0  fields  are  also  observable  in  the  sea  surface  data  set  and 
that  many  of  these  anomalies  can  be  more  strongly  associated  with  tectonic 
features  due  to  the  increased  resolution  of  the  sea  surface  field.  The  two 
data  sets  are  complementary,  in  areas  of  reduced  track  density  the  MAGSAT 
field  provides  higher  resolution  and  the  MAGSAT  data  set  also  records  the  all 
important  vector  field  data.  In  areas  of  high  track  density  the  sea  surface 
data  can  be  used  to  refine  the  MAGSAT  data  processing  techniques  and  also 
provide  further  constraints  on  source  models. 

Harrison  and  Carle  (1981)  recently  examined  the  spectral  content  of 
several  long  profiles  over  the  Pacific  and  Atlantic  oceans.  They  concluded 
that  intermediate  wavelength  energy  was  observable  in  marine  magnetics  pro¬ 
files.  Shure  and  Parker  (1981)  introduced  a  caveat  in  such  studies  of  long 
profiles.  Spectral  analyses  of  long  profiles  assume  that  all  anomalies  are 
lineated  perpendicular  to  the  track.  The  spectral  analyses  are  therefore 
highly  susceptible  to  three  dimensional  anomalies  or  oblique  lineations  which 
fold  energy  into  the  longer  wavelength  portions  of  the  spectrum.  For  conve¬ 
nience,  we  assign  the  term  'strike-aliasing'  to  this  phenomenon.  Shure  and 
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Parker  (1981)  concluded  their  study  by  bandpass  filtering  the  magnetic  anomaly 
pattern  in  an  area  of  the  Juan  de  Fuca  Ridge.  The  area  was  chosen  because  of 
its  dense  magnetic  surveys  perpendicular  to  the  strike  of  the  seafloor  anomaly 
pattern.  Shure  and  Parker  concluded  that  intermediate  wavelength  anomalies 
were  not  present  within  their  study  area.  The  MAGSAT  field  and  our  study 
(Figure  12)  both  shew  that  the  region  of  the  Juan  de  Fuca  ridge  is  a  region  of 
subdued  magnetic  anomalies  in  the  intermediate  wavelengths  which  strike  per¬ 
pendicular  to  the  track  direction.  Therefore,  it  cannot  be  concluded  from  the 
Shure  and  Parker  study  that  strike  aliasing  is  the  cause  of  the  observed 
intermediate  wavelength  field  in  sea  surface  profiles. 

Our  study  is  designed  to  circumvent  the  problem  of  'strike-aliasing'  by 
filtering  randomly  oriented  profiles  over  large  areas.  Our  technique  is  simi¬ 
lar  to  one  adopted  by  Nomura  (1979)  in  a  study  of  the  North  Pacific  interme¬ 
diate  wavelength  field  but  differs  in  methods  of  external  field  correction, 
data  density  and  areal  coverage. 

Figure  1  shows  the  results  of  this  study,  i.e.  the  intermediate  wave¬ 
length  total  field  anomalies  obtained  from  filtered  marine  anomaly  profiles 
over  the  North  Pacific.  Figure  2  displays  the  MAGSAT  anomaly  field  (Lange  1  et 
al.,  1982)  observed  at  400  km.  above  the  same  region  as  Figure  1.  Though 
there  is  some  difference  between  the  two  fields  of  Figs.  1  and  2  in  the  long 
wavelength  regional  fields,  close  examination  of  the  two  fields  reveals  a 
remarkable  correlation  which  will  be  discussed  later.  Figure  3  shows  the 
secular  variation  of  the  total  field  obtained  from  the  marine  data  set  for 
Epoch  1970.  The  intermediate  wavelength  anomaly  field  displays  a  strong 
correlation  to  bathymetric  features  (Fig.  4)  which  will  help  in  the  isolation 
of  source  bodies  and  the  development  of  improved  analysis  techniques  for 
MAGSAT  and  future  satellite  magnetic  surveys. 
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The  Data  Set: 

The  data  set  used  in  preparing  the  magnetic  anomaly  maps  (Fig.  1  and  Fig. 

3)  were  obtained  from  Lamont -Doherty  data  files  and  the  National  Geophysical 
and  Solar  Terrestrial  Data  Center  (NGSDC ) .  No  aeromagnetic  data  were  uti¬ 
lized.  The  data  set  consists  of  370  ships  tracks  and  over  38,000  filtered 
data  points  spanning  approximately  1,500,000  nautical  miles.  Figure  5  dis¬ 
plays  the  track  density  incorporated  in  this  study.  In  Figure  6  we  show  a 
histogram  of  the  data  coverage  with  respect  to  time.  The  data  coverage  spans 
the  years  1960  to  1980  with  a  mean  about  1970.  Therefore  we  will  define  our 
intermediate  wavelength  and  secular  variation  fields  to  epoch  1970.  Figure  6 
might  also  serve  as  an  indication  of  the  variation  in  funding  for  marine  geo¬ 
physical  exploration  during  the  past  two  decades. 

The  Technique: 

The  marine  data  measure  magnetic  fields  from  several  sources  of  internal 
and  external  origin.  Among  the  external  sources  are  magnetic  storms  and  the 
associated  DS  and  DST  variations,  Auroral  and  Equatorial  electrojets,  and 
ionospheric  SQ  currents  Nagata  and  Fukushima  (1971).  Other  sources  such  as 
magnetospheric  boundary  currents  or  micropulsations  are  either  of  low  ampli¬ 
tude  or  too  high  a  frequency  to  affect  our  study.  Internal  sources  include  the 
core  field  and  shallow  crustal  magnetic  distributions  of  structural  and  sea 
floor  spreading  type  caused  by  induced  or  remanent  magnetization.  Lithosphe 
ric  scale  magnetization  distributions  are  thought  to  be  the  sources  of  the 
intermediate  field.  Our  objective  in  filtering  the  surface  field  is  the  remo 
val  of  the  anomaly  fields  due  to  all  sources  except  those  of  lithospheric 
scale  and  intermediate  wavelength  within  the  interval  3000  to  300  km. 
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All  marine  data  are  first  reduced  to  anomaly  form  using  a  Pre-MAGSAT  ref¬ 
erence  field  of  degree  13  furnished  by  Dr.  R.  Langel  -MAGSAT  Project  Scien¬ 
tist.  The  model  includes  secular  variation  estimates  also  to  degree  13.  The 
field  model  serves  as  a  nominal  regional  field  which  can  easily  be  substituted 
later  in  the  analysis.  All  variations  in  the  potential  field  for  wavelengths 
greater  than  3000  km.  are  therefore  removed  from  the  observed  data  if  the 
PreMAGSAT  field  model  is  correct. 

The  total  field  anomalies  of  wavelength  less  than  300  km  were  filtered 
from  the  profiles  using  a  Gaussian  filter.  The  random  filtered  data  points 
were  then  processed  within  2x2  and  6x6  degree  bins.  The  random  orientation  of 
the  tracks  and  the  areal  nature  of  the  study  should  overcome  the  effects  of 
strike  aliasing.  A  true  test  of  the  effectiveness  of  our  technique  lies  in 
the  coherence  of  the  final  data  set  and  its  comparison  to  the  MAGSAT  anomaly 
field. 

We  have  applied  a  Gaussian  weighted  filter  to  all  anomaly  profiles  to 
eliminate  aliasing  from  shorter  wavelength  crustal  anomalies.  The  filter  is 
adapted  from  McKenzie  et  al.  (1980)  and  has  a  spectral  cutoff  of  300  km.  We 
have  filtered  point-by-point  values  of  the  magnetic  anomaly  M  (a-,tj),  compu¬ 
ted  Solar  Quiet  Variation  SQ(aj,tj),  and  DST  estimates  DST(ij,t-). 

where  Sj  =  The  position  vector  of  data  point 
tj  =  the  time  of  data  point 


The  Filtered  Magnetic  anomaly  becomes: 


A(b.,t.)  =  EM(a.,t.)  u>./E  to. 
i  J-  j  J  J  J  J 


(1) 
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where: 


b .  =  £  a .  m./Z  in . 

1  j  J  J  j  J 


2  2 

u).  =  exp(  — x  -  / cr  ) 
j  J 


A(&£,ti)  =  Averaged  anomaly  value  relocated  to  point 
(J  =  half  width  of  the  Gaussian  filter 

Xj  =  is  distance  to  point  aj  from  the  filter  centroid. 
The  spectral  response  of  the  Gaussian  antialiasisng  filter  is 


2  2 

/  n  \  \  _  (-S  O  /4) 

\2)  F(s;  =  exp 

where  s  =  scalar  wave  number. 

Figure  7  shows  the  application  of  the  filter  to  a  given  profile.  The 
filter  is  applied  by  moving  the  centroid  of  the  filter  at  a  given  increment 
along  the  ship's  track.  All  data  within  a  given  radius  of  the  centroid  is 
averaged  and  the  resultant  value  is  assigned  to  the  centroid.  The  filter  is 
incremented  at  50  kilometer  intervals  along  track  to  eliminate  aliasing.  Note 
that  the  filtered  value  is  assigned  to  a  location  defined  by  unit  vector  b^ 
which  is  not  necessarily  coincident  with  the  filter  centroid.  The  Gaussian 
half  width  of  a  =  100  km  was  chosen  to  give  a  low  pass  cutoff  of  approximately 
300  km  wavelength.  Figure  8  displays  the  spectral  response  of  the  Gaussian 
filter  applied  to  the  magnetic  data  set.  Figure  4  also  displays  the  expected 
amplitudes  and  wavelengths  of  various  internal  and  external  sources  to  be 
measured  by  a  shipboard  magnetometer  moving  at  10-15  knots.  Note  that  DST, 
DS,  SQ,  and  the  crustal  sources  are  of  sufficient  amplitude  to  cause  consider¬ 
able  error  if  these  field  components  are  not  properly  filtered. 
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The  marine  magnetic  data  are  scalar  total  field  measurements.  The  ano¬ 
maly  field  is  obtained  by  removing  a  reference  field  model  determined  from  a 
spherical  harmonic  expansion.  The  assumption  is  made  that  the  observed  field 
vector  is  in  the  direction  of  the  reference  field  vector.  This  can  be  ex¬ 
pressed  as: 

(3)  T(x,y)  =  T(x,y)  -  T(x,y) 

ANO  OBS  REF 

-  ~t (x, y)  •  V  A(x, y)  +  t(x,y)  •  V  A(x,y) 

REF  OBS  REF  REF 

(4)  T(x, y)  3  -t(x,y)  •  V  A(x,y) 

ANO  REF  ANO 

A(x,y)  =  Magnetic  scalar  potential 

T(x,y)  =  Magnitude  of  magnetic  field  vector 

t(x,y)  =  £(x, y)  i  +  ra(x,y)  j  +  n(x,y)k 

=  unit  vector  in  the  direction  of  the  magnetic  field  vector 

A  A 

tOBS  5  fcREF 

x,y,z  is  the  coordinate  system  pointing  in  the  North,  East, 
and  down  direction  respectively  with  direction  cosines  l,m,n. 
Subscripts:  OBS,  REF,  ANO  refer  to  the  observed,  reference  and  anomaly 
fields  respectively. 

Taking  the  two  dimensional  Fourier  Transform  of  equation  4  we  obtain 

(5)  T(u,v)  =  (  -1)  [n*s  A(u,v)  +  i(2*uA(u,v)  +  m*vA(u,v))] 

ANO  ANO  ANO  ANO 

where  s  =  (u^+v^)*^ 

=  scalar  wave  number 
*  is  the  convolution  operator 
Applying  the  anti-aliasing  filter  to  equation  3  we  obtain 

2  2 

-(c  s  ) 

F(u,v)  =  T(u,v)  e 
ANO 
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Note  that  if  the  magnetic  field  vector  is  not  constant  over  the  study 
area,  the  spectrum  of  the  total  field  anomaly  will  be  convolved  or  smoothed  by 
the  spectrum  of  the  total  field  direction  cosines.  The  effect  is  to  smear  the 
spectral  energy  to  adjacent  wavelengths.  For  small  study  areas  =106  km2, 
t(x,y)  can  be  thought  of  as  constant  and  our  filtering  process  is  a  bandpass 
for  3000  300  km.  Components  of  the  field  are  not  subject  to  this  convolu- 
tion  effect.  Therefore,  some  intermediate  wavelength  anomalies  may  appear  in 
the  total  field  pattern  and  not  in  the  component  field. 

Fig.  9  displays  the  correlation  between  DS,  DST,  and  Ap  (Suguira  and 
Chapman,  1960,  Mayaud,  1980) .  It  can  be  seen  that  during  a  large  magnetic 
storm,  DS  and  Ap  fall  off  at  approximately  equal  rates  after  the  main  phase  of 
the  storm.  This  may  be  due  to  the  relationship  between  substorms  and  magneto- 
spheric  convection.  The  compilations  of  Fig.  9  suggest  that  the  contributions 
of  DS  can  be  minimized  by  excluding  data  acquired  when  Ap>30.  The  remaining 
ring  current  field,  or  DST,  can  be  removed  using  the  estimates  of  Sugiura 
(1963)  available  from  the  NGSDC ,  when  Ap  <  30,  although  the  problem  of  annual 
variation  remains  unresolved  for  the  moment. 

Figure  10  is  a  histogram  of  the  occurrence  of  Ap  and  DST  values  for  our 
data  set.  Less  than  5%  of  the  utilized  data  fall  in  the  Ap  >30  range.  The 
effect  of  the  auroral  electrojet  was  avoided  by  restricting  the  analysis  to 
latitudes  less  than  50°N.  The  equatorial  electrojet  is  also  observable  in  the 
marine  data  within  the  region  ±5°  from  the  magnetic  equator  (Handschuraacher, 
1976).  We  have  not  restricted  our  data  in  time  or  area  to  avoid  the  effect  of 
the  equatorial  electrojet;  therefore  care  mist  be  taken  in  interpretations 
near  the  magnetic  equator.  The  location  of  the  magnetic  equator  is  shown  in 
as  a  stipled  line  in  Figures  1,2,3,14,  and  15. 
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The  diurnal  variation  or  SQ  can  provide  a  significant  level  of  variation 
in  the  observed  anomaly  pattern.  We  have  removed  SQ  utilizing  the  model  of 
Malin  (1973).  Figure  11  (a,b,c,d,e)  display  the  SQ  model  values  and  the 

variation  in  filtered  anomaly  data  with  respect  to  local  time.  A  total  of 
over  38,000  filtered  data  points  were  examined  and  the  number  of  data  points 
for  each  band  is  shown  in  the  figures.  The  diurnal  variation  is  clearly 
observed  in  the  data  set  and  on  the  average  the  Malin  model  effectively 
removes  this  field.  Yearly,  seasonal  and  daily  variations  in  SQ  were  not  con¬ 
sidered.  In  performing  the  averaging  it  was  noticed  that  the  mean  of  the  ano¬ 
maly  data  showed  a  latitudinal  dependence.  Figure  12  displays  the  plot  of 
filtered  anomaly  means  and  SQ  means  with  respect  to  latitude.  The  variation 
in  the  mean  anomaly  is  presumably  due  to  errors  in  the  low  order  terms  of  the 


PreMAGSAT  regional  field. 

The  ship's  field  at  the  sensor  is  usually  less  than  5-15  nT ,  Bullard  and 
Mason  (1963),  and  varies  with  the  ship's  heading.  Since  we  are  using  randomly 
oriented  tracks  within  a  given  study  area  we  conclude  that  the  effect  is  sub¬ 
dued  in  a  large  data  set.  A  more  serious  problem  is  the  effect  of  poor  navi¬ 
gation.  Inaccurate  navigation  results  in  the  calculation  of  an  incorrect  re¬ 
ference  field  which  may  be  several  tens  of  gammas  in  error.  Inspection  of 
individual  data  sets  and  the  culling  of  problematic  data  have  reduced  the 
effect  of  poor  quality  data  in  our  final  compilation. 

The  filtered  values  are  assembled  within  2  x  2  and  6  x  6  bins.  The 
sizes  of  the  bins  were  chosen  to  optimize  stability  of  the  solutions  and  spa¬ 
tial  resolution.  The  anomaly  values  within  a  given  bin  are  considered  to  be  a 
function  of  DST,  SQ,  secular  variation  and  internal  intermediate  wavelength 


sources. 
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The  data  within  a  given  bin  are  approximated  by  the  -function 
A(t )  =  Ayo  +  Vt  +  DST  +  SQ 

where 

A  =  observed  anomaly  (intermediate  wavelength  bandpass) 

Ajq  =  anomaly  at  epoch  1970 

V  =  secular  variation  anomaly  estimated  over  the  6°x6°  grid 
t  =  time  in  years  since  1970 

The  application  of  the  DST  and  SQ  estimates  significantly  improves  the 
regression  within  most  bins.  Figure  13  displays  a  bin  centered  on  23°  North, 
171°  East.  Note  the  improvement  in  the  distribution  of  filtered  anomaly  data 
with  respect  to  the  final  linear  regression  for  the  corrected  bin  values. 
Corrections  for  SQ  tend  to  reduce  the  dispersion  of  filtered  anomaly  values 
for  a  given  track  line:  the  DST  correction  generally  improves  the  regression. 

Total  Field  Secular  Variation  Anomalies: 

The  secular  variation  anomalies  (V)  were  computed  as  described  above  by 
linear  regression  within  6°  x  6°  bins.  A  secular  variation  anomaly  is  the 
error  in  estimate  of  secular  change  of  the  total  field  after  the  removal  of 
the  Pre-MAGSAT  model  (n  <  13).  Significant  secular  variation  anomalies  in  the 
total  field  are  seen  at  wavelengths  shorter  than  5000  km  and  amplitudes  as 
large  as  60  nT/yr  west  of  Hawaii.  (Figure  14)  .  These  anomalies  show  strong 
spatial  correlation  to  the  total  field  anomaly  pattern  in  Figure  15  which  is 
referenced  to  the  pre-MAGSAT  field.  For  example  compare  the  anomalies  cen¬ 
tered  on  30N,  175E  and  3 ON,  165W  in  Figures  14,  15  and  10.  This  can  easily  be 
observed  by  superimposing  the  two  maps. 

Alldredge  et  al.  (1963),  Harrison  and  Carle  (1981)  and  Carle  and  Harrison 
(1982)  have  shown  that  the  low  order  terms  in  the  potential  field  expansions 


may  contribute  to  high  order  terms  in  the  total  field.  Therefore,  the  ob¬ 
served  secular  variation  anomalies  may  reflect  errors  in  the  secular  variation 
estimates  of  the  Pre-MAGSAT  model  for  n<13.  Because  of  suspected  inaccuracies 
in  the  original  reference  field,  we  have  calculated  the  total  observed  fil¬ 
tered  field  and  its  associated  secular  variation  by  adding  the  field  estimates 
of  the  original  model  field  to  the  2°  x  2°  and  6°  x  6°  grid  values  for  total 
field  and  secular  variation  respectively.  The  anomaly  field  of  Figure  1  was 
calculated  by  removal  of  DGRF  1970  Field  (IAGA,1981)  to  order  10.  Figure  2  is 
a  plot  of  the  total  observed  secular  variation,  note  that  the  5000  km  wave- 
length  features  have  disappeared.  Some  localized  short  wavelength  features  do 
remain  in  the  secular  variation  field.  These  features  are  likely  artifacts  of 
strong  local  anomalies  and  insufficient  data  coverage. 

Core  phenomena  might  also  contribute  to  secular  variation  for  wavelengths 
shorter  than  3000km.  Lowes  (1974)  notes  that  the  energy  in  the  secular  varia¬ 
tion  appears  to  fall  off  less  rapidly  than  the  stationary  components.  This  is 
likely  due  to  the  westward  drift  of  the  nondipole  field  components.  There¬ 
fore,  the  remaining  short  wavelength  secular  variation  anomalies  may  be  due  to 
high  order  core  field  phenomena.  This  is  interesting  in  that  we  may  be  able  to 
identify  some  of  the  intermediate  wavelength  anomalies  as  high  order  core 
field  phenomena. 

Finally  we  must  consider  that  the  abrupt  truncation  of  a  reference  field 
at  a  given  order  will  generate  an  oscillation  or  ringing  (the  Gibbs  phenoraena- 
Chapman  and  Bartels,  p.  561,  1940;  Bracewell,  1971)  in  the  anomaly  field  for 
wavelengths  longer  than  the  high  pass  cut-off  in  this  case  \>  3000  km  or  n>13. 
Since  the  reference  field  for  both  Figures  15  and  2  (PreMAGSAT  reference 
field)  and  the  MAGSAT  and  Pogo  fields  were  truncated  at  N=13  we  would  expect 
similar  ringing  phenomena.  However,  the  DGRF  reference  field  was  truncated  at 
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order  and  degree  10  therefore  the  ringing  phenomena  should  differ  between 
these  fields.  Similarities  between  the  DGRF  field  and  the  MAGSAT  field  should 
therefore  have  physical  significance  since  we  can  eliminate  the  Gibbs  effect 
as  a  cause  for  anomaly  correlations.  The  bandpass  for  figure  1  after  aver¬ 
aging  with  2°x2°  bins  and  removal  of  the  degree  10  DGRF  is  4000  km  to  400  km. 

Intermediate  Wavelength  Total  Field  Anomalies: 

In  general  a  strong  correlation  exists  between  the  MAGSAT  anomaly  field 
Figure  2  and  the  filtered  marine  data  Figure  1  for  the  North  Pacific.  Because 
the  two  data  sets  are  derived  from  independent  data  sets  and  different  reduc¬ 
tion  techniques  were  utilized  for  each  field.  The  correlation  therefore 
suggests  that  the  techniques  applied  to  each  data  set  are  appropriate  and  that 
the  observed  anomalies  are  not  an  artifact  of  processing.  Positive  anomalies 
are  observed  in  the  satellite  and  surface  fields  along  the  Western  Pacific 
Trench  system,  (35N,  130E),  the  Indonesian  islands  (5N,  130E)  and  the  Central 
Pacific  (25N,  170E).  The  northwestern  Pacific  is  characterized  by  high 
amplitude  intermediate  wavelength  anomalies  while  the  northeastern  Pacific 
exhibits  a  more  subdued  pattern.  The  basement  relief  in  the  North  Pacific  is 
generally  greater  in  the  northwestern  Pacific  than  the  northeastern  Pacific 
(Figure  4) .  The  basement  relief  northwestern  Pacific  is  dominated  by  lineated 
seamount  chains  attributed  to  hot  spot  trends  or  fracture  zones  (Wuksibm  1965; 
Morgan,  1973) . 

A  first  look  at  the  two  intermediate  wavelength  fields  shows  that  the 
horizontal  gradients  of  the  sea  surface  data  are  generally  steeper  than  those 
of  the  MAGSAT  field.  This  is  expected  since  the  MAGSAT  field  should  be  equi¬ 
valent  to  the  upward  continued  sea  surface  data  provided  no  errors  exist  in 
either  data  set.  The  gradients  of  the  sea  surface  data  generally  outline  tec- 
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tonic  features  such  as  the  western  Pacific  trench  system,  the  mid  Pacific 
Seamounts,  the  Hess  and  Shatsky  Rises,  the  margins  of  the  western  Pacific  Back 
Arcs  Basins,  the  Hawaiian  Emperor  Seamount  Chain,  and  to  a  lesser  extent,  the 
Eastern  Pacific  Fracture  Zones  such  as  the  Clipperton,  the  Clarion  and  the 
Mendocino. 

Broad  anomalies  ranging  from  1000  km  to  greater  wavelengths  are  also  ob¬ 
served  in  both  fields.  The  most  prominent  of  these  we  chose  to  call  the 
Emperor  Anomaly.  The  Emperor  Anomaly  at  14  nT.  ,  is  the  largest  amplitude 
anomaly  observed  in  the  MAGSAT  field  over  the  ocean  basins.  Our  filtered  sea 
surface  data  set  also  show  this  anomaly  with  a  maximum  amplitude  of  150  nT. 
focused  at  the  bend  of  the  Emperor  Hawaiian  seamount  chain  (Figures  1,2,  and 

16  to  19)  where  it  intersects  the  Hess  Rise. 

The  Line  Island  Chain  and  Mid-Pacific  seamounts  separate  the  broad  posi¬ 
tive  regional  anomaly  (the  Emperor  Anomaly)  to  the  north  from  a  generally 
negative  residual  field  to  the  south  (Figure  1).  The  anomaly  field  to  the 
east  of  the  Pacific  trench  system  is  generally  negative  as  are  the  Miocene 
back  arc  basins  from  the  Sea  of  Japan  to  the  Parece  Vela  Basin.  The  older 
Eocene  age  Philippine  Basin  displays  a  relative  positive  anomaly  field 
(Figures  1,  2  and  23). 

The  magnetic  pattern  over  the  Emperor  anomaly  of  the  Western  Pacific  is 
observed  in  both  the  sea  surface  data  and  the  satellite  field.  Figures  16  and 

17  display  the  pattern  of  the  intermediate  anomalies  with  respect  to  regional 
bathymetry.  The  Emperor  is  bounded  on  its  southern  and  western  margins  by  a 
steep  negative  gradient  of  20-40YnT/100km.  The  southern  boundary  is  strongly 
correlated  to  the  Mid  Pacific  Seamount,  the  Western  Boundary  also  has  a  sea¬ 
mount  lineament  and  the  Shatsky  Rise  to  the  Northwest.  The  Northeastern  Boun¬ 
dary  is  likely  the  Hess  Rise  though  the  maximum  of  the  Emperor  Anomaly  is 
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located  above  the  bend  in  the  Emperor  Seamount  Chain,  and  the  intersection  of 
the  Mendocino  Fracture  Zone. 

Though  the  sharper  gradients  appear  to  be  due  to  morphologic  features, 
the  broad  positive  magnetic  anomaly  of  the  Emperor  Anomaly  appears  to  be  rela¬ 
ted  to  a  zone  of  high  amplitude  seafloor  spreading  type  anomalies  of  Mesozoic 
age.  In  Figures  18  and  19  we  show  the  MAGSAT  and  filtered  seasurface  fields  in 
relation  to  the  seafloor  spreading  anomalies.  Note  that  the  maximum  in  the 
intermediate  fields  appears  to  correlate  strongly  with  the  maximum  amplitude 
of  the  seafloor  spreading  anomalies.  The  amplitude  of  the  seafloor  spreading 
anomalies  may  reflect  an  enhancement  in  the  magnetic  mineralogy  of  the  litho¬ 
spheric  plate. This  suggests  to  us  that  in  addition  to  the  magnetic  fields  of 
the  seamount  chains,  the  Hess  and  the  Shatsky  Rises,  the  region  of  the  Emperor 
Anomaly  may  also  be  associated  with  a  lithosphere  of  generally  hgher  suscepti¬ 
bility  . 

Overlaying  the  magnetic  anomaly  maps  1  and  2  over  the  basin  bathymetric 
shows  that  linear  negative  anomalies  (40-50  nT  negative)  are  located  above  the 
older  seamounts  of  the  region  (Figures  16  and  17) .  The  seamounts  of  the  Mid 
Pacific  are  generally  Mid  to  Late  Cretaceous  in  age  and  were  formed  during  the 
Cretaceous  normal  period  south  of  the  magnetic  equator  (Thiede  et  al.  ,  1981, 
Lancelot,  1978;  Lancelot  and  Larson,  1975).  The  Emperor-Hawaiian  chain  are 
Late  Cretaceous  to  Cenozoic  in  age  and  were  formed  near  20° N  (Kono,  1980, 
Harrison  et  al. ,  1975).  The  northward  motion  of  the  Pacific  has  moved  them  to 
their  present  location  (Fig.  20). 

Figures  16  and  17  show  that  the  aggregate  magnetic  anomaly  pattern  of 
these  seamounts  is  strongly  negative  at  sealevel  and  weakly  negative  at  the 
satellite  altitudes.  A  comparison  of  Figures  16  and  17  shows  the  strong 
correlation  between  the  MAGSAT  and  sea  surface  fields.  The  figures  also  show 
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that  decreased  sensitivity  of  MAGSAT  to  north-south  striking  anomalies  due  to 
its  orbital  inclination.  Note  that  the  Emperor  Seamount  anomaly  is  not  ob¬ 
served  in  MAGSAT  while  the  Emperor  Trough  anomaly  is  apparent  in  both  fields. 

Figures  16  and  17  shew  that  a  positive  anomaly  exists  of  the  Northwest 
and  Northeast  corners  of  the  Emperor  anomaly.  These  appear  to  be  associated 
with  zones  of  thickened  crust  named  the  Shatsky  and  Hess  Rises. 

Because  most  large  seafloor  features  X>  100  km  are  isostatically  compen¬ 
sated,  we  expect  that  intermediate  wavelength  topography  will  be  reflected  in 
mantle  topography.  Wasilewski  et  al.  (1979)  have  suggested  that  the  Moho 
forms  a  lower  magnetic  boundary  within  the  lithosphere.  The  work  of 
Wasilewski  et  al.  is  based  on  susceptibility  measurements  of  ultramafics.  A 
direct  relationship  should  therefore  exist  between  seafloor  topography  and 
intermediate  wavelength  anomalies  if  the  hypothesis  of  Wasilewski  et  al. 
(1979)  is  correct  and  if  oceanic  relief  is  isostatically  compensated  by  litho¬ 
spheric  flexure  and  crustal  thickening  (Watts,  1978). 

The  question  of  whether  the  observed  fields  are  due  to  induced  or  rema¬ 
nent  magnetization  is  especially  interesting  since  an  induced  anomaly  of  the 
amplitude  observed  over  the  Mid  Pacific  Seamounts  would  require  a  substantial 
root  and  would  support  Wasilewski's  hypothesis.  On  the  other  hand  a  remanent 
magnetic  source  could  significantly  aid  in  determining  absolute  plate  rota¬ 
tions  from  paleoraagnetic  data.  In  the  region  covered  by  Figures  16,17  we  have 
a  good  record  of  plate  motion  and  paleomagnet ics  as  recovered  by  the  Deep  Sea 
Drilling  Program  and  the  work  of  such  authors  as  Harrison  et  al. ,  1975;  Kono, 
1980;  Lancelot  and  Larson,  1975;  Lancelot,  1978;  Thiede  et  al.  ,  1981;  and 
Vallier  et  al.  ,  1981.  In  order  to  determine  the  most  important  component  of 
magnetization  we  have  formulated  two  simple  models.  The  first  model  is  a  2- 
dimensional  prism  of  thickness  4  km  and  magnetization  3.3  A/m  (.0033  entu/cc) 
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and  polarity  according  to  references  previously  mentioned  and  listed  in  the 
figure  caption  of  Figure  21.  The  second  model  is  an  inductively  magnetized 
prism  of  thickness  13  km  and  susceptibility  4xl0“2  SI  (.0033  emu/cc).  The 
upper  surface  of  each  prism  is  at  2  km  depth.  Figure  21  displays  the  observed 
anomaly  at  seasurface  and  at  MAGSAT  altitudes.  A  comparison  of  the  model  pro¬ 
files  to  the  various  seamount  chains  and  rises  of  figures  16  and  17  shows  that 
only  for  the  Marshall-Gilbert,  Mid-Pacific  and  Emperor  Seamounts  is  it  possi¬ 
ble  to  distinguish  between  the  two  models.  The  induced  anomaly  model  for  the 
Marshall-Gilbert  and  Mid  Pacific  seamounts  is  much  more  highly  skewed  than  the 
observed  data.  The  observed  data  match  the  remanent  model  quite  well  with  the 
with  the  minimum  of  the  anomaly  directly  over  the  center  of  mass.  The  Emperor 
Seamounts  strongly  support  a  remanent  model  since  the  induced  model  requires  a 
positive  anomaly  where  a  negative  (-10  y)  anomaly  is  observed  in  the  filtered 
seasurface  field.  Unfortunately  the  North-South  lineated  anomaly  is  only 
slightly  recorded  in  the  MAGSAT  field.  This  is  presumably  due  to  the  direc¬ 
tional  sensitivity  of  the  MAGSAT  field. 

Figure  23  displays  the  filtered  sea  surface  field  with  respect  to  the 
tectonic  features  of  the  Western  Pacific  island  arc  system.  Contours  are 
extended  over  islands  and  onto  land  areas.  This  is  a  result  of  our  numerical 
processing  and  represents  an  interpolation  or  extrapolation  of  the  marine  data 
set.  Therefore  care  mist  be  exercised  in  interpreting  these  results.  The 
correlations  between  the  filtered  sea  surface  data  and  the  tectonic  structure 
of  the  region  are  surprisingly  strong.  The  similarity  of  the  filtered  sea 
surface  field  to  the  MAGSAT  field.  Figure  24,  is  generally  poor.  Some  simi¬ 
larity  exists  along  the  Japan  trench  and  over  the  island  arcs  but  the  corre¬ 
lation  worsens  over  the  back  arc  system.  This  may  be  due  to  the  insensitivity 


of  the  MAGSAT  field  to  north-south  lineations  which  pre-dominate  in  the 
region. 

One  can  see  a  strong  correlation  in  Figure  23  to  virtually  all  known 
tectonic  lineaments.  These  correlations  include  local  positive  anomalies 
bordering  all  trenches  including  the  Kuril,  Japan,  Bonin,  Ryukyu,  and  the 
Philippine  trenches.  The  only  exception  being  the  Mariana  trench.  A  local 
negative  follows  the  abandoned  spreading  axis  of  the  Japan  Basin,  Shikoku  and 
Parece  Vela  Basins.  These  basins  are  likely  Miocene  to  Oligocene  in  age 
(Mr ozowski  and  Hayes,  1979;  Karig,  1971,  1975;  Weissel,  1981;  Karig  et  al.  , 
1975).  The  older  (Eocene)  West  Philippine  Basin  (Mrozowzki  et  al. ,  1982) 

shows  a  local  positive  anomaly.  A  weak  magnetic  gradient  is  observed  over  the 
Central  Basin  Fault  and  stronger  gradients  are  located  along  the  boundaries  of 
the  Philippine  Basin. 

Though  it  is  quite  obvious  that  the  intermediate  wavelength  anomalies 
appear  to  separate  regions  of  the  Western  Pacific  with  differing  tectonic 
histories  the  causes  of  the  magnetic  contrast  are  not  obvious  and  the  appli¬ 
cation  of  a  single  model  is  fraught  with  contradictions. 

The  general  correlation  of  moderate  heat  flow  to  the  Japan-Parece  Vela 
Basin  negative  anomaly  (Yanagisawa  et  al. ,  1982)  is  contradicted  by  the  older 
(Eocene)  and  lower  heat  flow  of  the  Oki-Daito  region.  Furthermore,  maximum 
heat  flow  is  observed  near  the  island  arc  system  and  not  near  the  axis  of  the 
negative  magnetic  anomaly  in  the  Parece  Vela  Basin  (Anderson  et  al.  ,  1978). 
Figure  25a  displays  a  plot  of  heat  flow  values  with  respect  to  magnetic  ano¬ 
maly.  Heat  flow  values  were  taken  from  Anderson  et  al. ,  1978.  Error  bars  on 
the  heat  flow  values  are  1  standard  deviation.  Of  course  the  heat  flow  values 
may  not  accurately  reflect  the  true  temperature  of  the  lithosphere  due  to  com- 
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plications  such  as  local  hydrothermal  circulation  or  topographic  effects.  The 
correlation  is  nonetheless  weak. 

Variations  in  the  chemistry  of  the  oceanic  lithosphere  may  reflect 
changes  in  the  bulk  magnetic  mineralogy  of  the  basins.  However,  Deep  Sea 

Drilling  results  shew  that  the  composition  of  West  Philippine,  Shikoku,  and 
Parece-Vela  Basin  basalts  closely  resemble  the  mid-ocean  ridge  basalts  with 
only  a  slight  increase  in  alkalinity  from  the  Philippine  to  the  Parece  Vela 
basin  (Zakariadze,  et  al.,  1980;  Dick,  et  al.,  1980;  Marsh,  et  al.  ,  1980;  and 
Wood  et  al.  ,  1980).  For  the  moment  however,  the  data  are  too  sparse  and  the 
correlation  are  too  subtle  to  draw  any  conclusions  about  the  relationship 
between  the  magnetic  anomaly  pattern  and  the  chemical  composition  of  the 

oceanic  lithosphere. 

There  appears  to  be  moderate  correlation  between  crustal  thicknesses  and 
magnetic  anomaly  amplitude  along  the  West  Mariana  Ridge  and  within  the  Parece- 
Vela  Basin  and  Japan  Sea  (Fig.  25b).  This  would  support  the  alternate  hypo¬ 
thesis  of  Yanagisawa,  et  al. ,  1982.  No  similar  correlation  is  observed  within 
the  West  Philippine  Basin  though  the  seismic  data  are  sparse  and  very  few  re¬ 
fraction  profiles  have  measured  the  depth  to  Moho  along  the  margine  of  the 

West  Philippine  Basin.  We  note  that  ultra  basic  have  been  dredged  from  within 
the  Parece  Vela  Basin  which  further  suggests  a  thin  crust  and  shallow  Moho 
(Mrozowski  and  Hayes,  1979). 

There  is  a  strong  suggestion  that  intermediate  wavelength  anomalies  deli¬ 
neate  tectonic  boundaries.  The  anomalies  are  moderately  correlated  to  crustal 
thickness,  and  weakly  correlated  to  the  measured  heat  flow  though  problems  in 
obtaining  proper  measurement  environments  may  have  biased  the  measured  heat 
flow  toward  lower  values  and  confused  the  correlation. 
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In  figure  26  we  show  that  a  satisfactory  model  for  the  western  Pacific 
can  be  formulated  by  varying  the  lower  boundary  of  the  magnetic  body  in  accor¬ 
dance  with  Moho  depth  as  derived  from  seismic  data  (Hayes  et  al.  ,  1978). 
Thermal  conditions  within  the  plate  may  further  alter  this  crustal  thickness 
mode  1. 

Lineated  negative  anomalies  associated  with  the  Mid  Pacific  seamounts  are 
observed  intersecting  the  Western  Pacific  trench  system  (Figures  1  and  23). 
It  is  curious  that  this  lineation  of  the  negative  anomalies  appears  to  conti¬ 
nue  into  the  back-arc  region.  It  is  interesting  to  speculate  that  the  subduc- 
tion  of  the  Seamount  chain  may  have  left  a  chemical  and  structural  imprint  in 
the  developing  back-arc  as  suggested  by  Kellerher  and  McCann  (1976). 

In  the  eastern  Pacific  there  is  a  remarkable  correlation  between  the 
major  fracture  zones  and  the  intermediate  wavelength  anomalies  observed  in  the 
MAGSAT  field  (Figure  2)  and,  to  a  lesser  extent,  in  the  2<,x2°  sea  surface 
field  (Figure  1).  Since  there  are  relatively  few  tracks  in  the  eastern 
Pacific  to  constrain  the  sea  surface  field,  the  MAGSAT  field  may  be  more  accu¬ 
rate  in  depicting  this  correlation,  particularly  for  east— west  lineated  anoma¬ 
lies.  In  Figure  27  we  show  the  geophysical  data  collected  by  the  C.V.  Hudson 
along  a  track  that  ran  up  longitude  150°W,  crossing  the  Clarion,  Molokai, 
Murray  and  Mendocino  fracture  zones.  (This  is  the  same  profile  studied  by 
Harrison  and  Carle,  1981.)  Large  magnetic  anomalies  are  observed  in  the  ob¬ 
served  and  smoothed  profiles  above  the  fracture  zones.  More  importantly,  it 
is  apparent  that  the  fracture  zones  correspond  to  changes  in  the  gradients  of 
the  2  x2  sea  surface  field  and  in  the  MAGSAT  field  as  shown  in  the  profiles 
in  the  upper  part  of  Figure  27. 

We  have  modeled  these  anomalies  as  due  simply  to  a  small,  .038  SI  (.003 
emu/cc),  susceptibility  contrast  across  the  fracture  zones  (Figure  28).  (The 
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anomalies  could  also  be  attributed  to  a  small  contrast  in  the  remanent  magne¬ 
tization  across  the  fracture  zones.)  The  fact  that  the  susceptibility  (or 
remanent  magnetization)  contrasts  are  maintained  over  long  distances  suggests 
that  they  may  represent  long  period  (30-50  m.y.)  variations  in  the  magnetic 
properties  of  the  crust  across  the  fracture  zones.  Long  term  variations  in 
basalt  chemistry  across  fracture  zones  has  been  observed  in  rock  samples 
collected  by  Challenger  drilling  in  the  Atlantic  Ocean  (Scientific  Staff, 
1982).  The  origin  of  these  long  term  variations  are  unknown,  but  may  reflect 
the  blocking  effect  of  transform  faults  on  longitudinal  flow  along  the  mid¬ 
ocean  ridge  system  or  as  a  manifestation  of  longitudinal  cellular  convection 
(Richter  and  Parsons,  1975). 

Conclusion: 

We  have  shown  that  in  regions  of  relatively  high  track  density,  it  is 
possible  to  recover  the  intermediate  wavelength  magnetic  anomaly  field  from 
marine  surveys.  Except  in  cases  of  north-south  lineated  anomalies  the  sea 
surface  anomaly  field  strongly  correlates  to  the  MAGSAT  intermediate  wave¬ 
length  field  over  marine  regions  therefore  verifying  the  validity  of  both 
fields. 

In  comparing  the  marine  data  to  the  MAGSAT  anomaly  pattern,  the  severity 

of  the  directional  sensitivity  has  become  apparent.  Though  we  cannot  quantify 

this  directionality  with  great  precision,  a  comparison  of  the  Emperor  seamount 

anomaly  in  the  sea  surface  and  the  MAGSAT  fields  suggests  that  the  attenuation 

of  north-south  anomalies  is  approximately  a  factor  of  4.  A  more  careful  spec¬ 
tral  analysis  of  the  two  fields  will  be  required  to  more  accurately  quantify 


this  effect. 
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The  observed  anomalies  over  the  central  correlate  strongly  with  tectonic 
lineaments  of  the  seafloor.  We  are  able  to  show  that  the  mid  Pacific  seamount 
province  and  rises  are  observable  in  the  satellite  field  and  that  at  least 
some  of  these  observed  anomalies  are  due  to  remanent  magnetization. 

The  western  Pacific  show  a  strong  correlation  between  tectonic  features 
of  the  island  arc  systems  and  observed  positive  anomalies.  This  correlation 
appears  to  be  due  more  to  variations  in  crustal  thickness  than  heat  flow.  The 
Pacific  plate  shows  many  intermediate  wavelength  magnetic  lineaments  which  may 
be  related  to  the  tectonic  history  of  the  Pacific  plate.  A  careful  analysis 
of  both  the  sea  surface,  and  MAGSAT  magnetic  fields  in  conjunction  with  the 
region's  geologic  data  will  undoubtedly  reveal  more  about  the  sources  of  these 

magnetic  anomalies  and  geologic  development  of  the  Pacific  and  other  oceanic 
plates. 
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1:  Filtered  seasurface  data  for  the  bandpass  4000  km  to  400  km  as  des¬ 
cribed  in  the  text.  The  anomaly  field  is  with  respect  to  the 
DGRF  1970  reference  field  (IAGA,  1981).  The  magnetic  equator  is 
indicated  as  a  stipled  line. 

2:  MAGSAT  anomaly  field  from  2°x2°  averages.  Average  altitude  is  400 
km  (Lange  et  al. ,  1982) . 

3:  Map  of  total  field  secular  variation  as  derived  from  the  filtered 
sea  surface  data.  The  magnetic  equator  is  indicated  by  a  stipled 
line. 

4:  Bathymetric  map  of  the  region  from  Chase  et  al.  ,  with  DSDP  drill 
sites  indicated. 

5:  Track  locations  of  the  data  set  used  in  this  study. 

6:  A  histogram  of  the  data  set  distribution  in  time.  Note  that  the 
data  set  is  normally  distributed  in  time  about  1970. 

7:  Application  of  the  Gaussian  filter  to  a  measured  magnetic  anomaly 
profile. 

8:  Spectral  Response  of  the  Gaussian  antialiasing  filter.  Marine 
crustal  sources  are  assumed  to  be  largely  due  to  seafloor 
spreading  type  anomalies  of  short  wavelength  and  high  amplitude. 

9:  A  correlation  between  the  temporal  variations  in  DS,  DST  and  Ap. 

10: Frequency  of  DST  and  Ap  values  in  the  marine  data  set. 

1 1 : Correlation  between  observed  diurnal  variation  and  the  SQ  model 
(Malin,  1973)  utilized  in  this  study.  Data  are  divided  into  10 


degree  latitudinal  bands. 
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Figure  12: Illustration  showing  the  the  PreMAGSAT  referenced  anomaly  field  of 
Figure  15  displays  a  strong  latitudinal  variation  presumably  due 
to  field  error.  The  SQ  model  is  stable  throughout  the  region. 
Diurnal  is  the  mean  anomaly  value  for  each  latitudinal  band  of 
Fig.  11.  SQ  is  the  mean  model  value  for  each  band. 

Figure  13:  Application  of  the  regression  analysis  to  data  within  a  2°x2°  km  at 
22°N  and  170°  East.  Note  that  corrections  for  DST  and  SQ 
significantly  improve  the  regression.  It  was  found  that  better 
stability  in  the  secular  variation  estimation  could  be  obtained 
for  6x6  degree  bins. 

Figure  14: Secular  variation  anomaly  field  estimated  from  6x6  degree  bins 
referenced  to  the  Premag  I  field.  Note  the  oscillation  in  values 
about  the  Hawaii  suggesting  an  error  in  the  model. 

Figure  15:6°x6°  anomaly  field  derived  from  the  PreMAGSAT  field.  Note  that  the 
major  secular  variation  anomalies  of  Figure  14  also  appear  in  the 
anomaly  field,  particularly  in  the  region  of  Hawaii.  This 
indicates  the  sensitivity  of  the  anomaly  field  to  the  selected 
reference  field. 

Figure  16:MAGSAT  field  contoured  at  a  .5nT  contour  interval.  Units  are 
.InT.  Note  that  there  is  a  slight  negative  above  the  Emperor 
Seamount  chain. 

Figure  17: Sea  surface  filtered  magnetic  anomaly  referenced  to  DGRF  1970 
field.  Contour  interval  is  10  nT.  Stipled  zones  are  regions 

above  the  4000  meter  isobath.  Note  the  strong  correlation  to  the 
MAGSAT  field  and  Figure  16.  Lineated  negative  anomalies  follow 
the  mid-Pacific  seamounts,  the  Marshall-Gilbert  and  Emperor  Sea- 
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mounts.  Positive  anomalies  overlie  the  Shatsky  and  Hess  Rises 
and  the  Emperor  Seamounts. 

Figure  18: MAGSAT  field  referenced  to  the  observed  seafloor  spreading  anomalies 
(Larson  and  Hilde,  1975). 

Figure  19: Filtered  sea  surface  field  with  respect  to  the  seafloor  spreading 
anomalies  (Larson  and  Hilde,  1975).  Note  that  the  seafloor 
spreading  anomalies  of  a  given  age  fall  off  in  amplitude  as  the 
Emperor  anomaly  decreases  in  amplitude.  We  take  this  as  evidence 
for  a  regional  variation  in  crustal  lithology  emplaced  during  the 
creation  of  this  portion  of  the  mid-Pacific  in  the  Jurassic  and 
Early  Cretaceous. 

Figure  20:?late  motion  diagram  for  DSDP  sites  as  calculated  by  Lancelot  and 
Larson  1975,  and  Lancelot  (1978). 

Figure  21:Model  calculations  for  uniformly  magnetized  prisms.  Two  models  are 
considered:  (1)  Induced  prism  200  km  wide  by  13  km  deep.  Sus- 

cept ibility=  .04  SI  (.0033  emu/cc);  (2)  Remanent  Model=  200  km 
wide  by  4  km  deep.  Magnet  ization=  3/3  A/m  (.0033  emu/cc).  The 
upper  of  each  model  is  at  2  km  depth.  Note  that  only  the  rema¬ 
nent  and  induced  models  of  the  Marshall  Gilbert/Mid— Pacific  and 
Emperor  Seamounts  are  significantly  different. 
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TABLE 


Azimuths,  remanent  and  present  field  vectors  are  as  follows: 


REMANENT  FIELD  PRESENT  FIELD 
STRUCTURE  AZIMUTH  IR  DR  Ip  Dp 


Emperor  Srats 

90° 

-42° 

0° 

52° 

3° 

Hawaiian  Smts 

30° 

31° 

0° 

o 

O 

7° 

Mid-Pacific  Smts 

0° 

-31° 

0° 

31° 

9° 

Marshall  Gilbert  Smts 

30° 

-31° 

0° 

31° 

9° 

Shatsky  Rise 

315° 

-10° 

0° 

e 

CM 

-2° 

Hess  Rise 

355° 

-10° 

0° 

48° 

8° 

Figure  22:  Conrad  1007  profile 

over  the 

Mid-Pacific 

seamounts  showing  upward 

cont inued 

field  and 

MAGSAT 

anomalies . 

and  the 

negative  anoraa ly 

associated  with  the  Mid-Pacific  Seamounts. 

Figure  23:  Filtered  sea  surface  anomaly  field  over  the  Western  Pacific,  45 

nT  have  been  added  to  the  observed  anomalies.  Note  the  strong 
correlation  between  the  region's  ridge  systems  and  the  magnetic 
field.  A  negative  follows  the  abandoned  back  arc  systems  of  the 
Japan  Basin,  Shikoku  Basin  and  the  Parece  Vela  Basin.  In  general 
a  positive  anomaly  is  observed  over  the  abandoned  and  present 
island  arc  systems.  The  only  exception  being  the  Mariana  Arc 
system. 


Figure  24:  The  MAGSAT  field  above  the  Western  Pacific.  The  correlation  is 

much  poorer  than  for  the  Central  Pacific.  Presumably  this  is  due 
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Co  the  directional  sensitivity  of  the  MAGSAT  field.  A  linear 
positive  is  observed  above  the  Arc  trench  system  but  the  detail 
of  the  sea  surface  field  is  not  observed. 

Figure  25a:  Plot  of  observed  heat  flow  averages  versus  magnetic  anomaly  com¬ 
piled  for  all  stations  from  Anderson  et  al.,  1970.  Heat  flow 
values  shew  a  broad  scatter  which  can  be  reduced  somewhat  by 
better  station  election.  However,  hydro  thermal  circulation  may 
have  seriously  biased  the  values  circulation  lower  heat  flow 
values . 

Figure  25b:  Correlation  between  magnetic  anomaly  value  and  the  depth  to  base¬ 
ment  as  taken  from  the  sonobuoy  compilations  of  two  correlations 
are  shown  one  for  the  Parece  Vela  Bas in-Sh ikoku-Japan  Basin  and 
the  other  for  the  West  Phillipine  Basin.  It  should  be  noted  that 
many  variables  will  affect  this  correlation  including  body  struc¬ 
ture  and  magnetization  distributions.  The  Parece  Vela  anomalies 
are  generally  linear  and  North-South  trending. 

Figure  26:  An  induced  model  based  on  the  compilations  of  the  western  Pacific 

seismic  data  may  explain  the  observed  intermediate  field  in  terras 
of  variations  in  the  Moho  depth.  The  susceptibility  contrast  of 
the  body  is  .063  SI  (.005  emu/cc).  The  subducting  plate  has 
little  the  effect  on  the  anomaly  field.  The  model  appears  dis¬ 
torted  is  horizontal  scale  due  to  large  vertical  exaggeration. 

Figure  27:  (Lower)  Magnetic,  gravity  and  bathymetry  from  the  C.V.  Hudson 

along  longitude  150°W.  (Upper)  Filtered  20x2°  sea  surface  data, 
upward  continued  to  400  km  (dashed),  and  the  MAGSAT  field  along 
the  same  track  as  the  Hudson  profile. 
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Figure  28:  Model  for  an  induced  anomaly  due  Co  a  susceptibility  contrast 
across  the  eastern  Pacific  fracture  zones.  Stippled  bodies  have 
a  .038  SI  (.003  emu/cc)  susceptibility,  white  bodies  have  a  zero 
susceptibility,  azimuth=0°,  Ip=50° ,  Dp=10° ,  field  strength=40000 
nT.  Filtered  profile  has  a  bandpass  between  400  and  2000  km. 


Appendix  A 


Tracks  and  filter  parameters  used  in  the  calculation  of  the 
sea  surface  intermediate  wavelength  anomaly  field.  Parameters 
and  track  ID's  refer  to  LDGO  program  and  data  libraries. 


APPENDIX  B 

Subroutine  FIELDG  and  the  associated  Pre-Magsat  coefficients 
used  in  the  generation  of  the  Pre-Magsat  reference  Field. 
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SUBROUTINE  FIELDG  (DLAT,  DLONG.  ALT,  TM,  NMX,  L,  X,  Y,  Z,  F) 

C  MODIFIED  JULY  22,  19S1 

FOR  DOCUMENTATION  OF  THIS  SUBROUTINE  AND  SUBROUTINE  FIELD  SEE  : 
NATIONAL  SPACE  SCIENCE  DATA  CENTER'S  PUBLICATION 

COMPUTATION  OF  THE  MAIN  GEOMAGNETIC  FIELD  FROM 
SPHERICAL  HARMONIC  EXPANSIONS 
DATA  USERS'  NOTE, NSSDC  68-11,  MAY  1968 
GODDARD  SPACE  FLIGHT  CENTER- GREENBELT, MD. 

EQUIVALENCE  (SHMIT( 1, 1 > , TG< 1, 1  )  > 

COMMON  NDPY  <  3  > , LPYR ( 4 ) , JBSYR 

COMMON  /CQEFFS/TG ( 18, 18), G(1S, 18), GT(18, 18), GTT(1B, 18), TLAST 
1, TZERO 

COMMON  /FLDCCM/ST , CT, SPH, CPH, R, NMAX , BT,  BP,  SR,  B 
DIMENSION  AID( 11 ), SHMIT( 18, 18) 

DATA  A/Q.  / 

A=6378.  16 

FLAT=1 .  -1.  /298.  25 

A2=A**2 

A4=A**4 

B2=(A*FLAT>**2 
A2B2=A2* ( 1 . -FLAT** 2 ) 

A4B4=A4*(1. -FLAT**4) 

c ******* IF  L.  LE.  0  DON'T  READ  COEFFICIENTS  FILE 

C*******IFLDATE1il  THE  SAME  AS  LAST  DON'T  UPDATE  COEFFICIENTS 
1  IF  (TM-TLAST)  17,19,17 


C*******READ 


FIELD 


1  =  1,  11 ) 


(213,  5X6HEP0CH=, F7.  1,  5X10A6, A3) 


74 


COEFFICIENTS 

READ  (3,3)  J,  K,  TZERO,  (AID(  I  ), 

FORMAT  (211,  IX,  F6.  1,  10A6,  A3) 

L=0 

FORMAT 
MAXN=0 
TEMP=0 

READ  (3,6,  END=74)  N,  M,  GNM,  HNM,  GTNM,  HTNM,  GTTNM,  HTTNM 
WR ITE ( 5, 6 ) N, M, GNM, HNM, GTNM, HTNM, GTTNM, HTTNM 
FORMAT  (213, 6F11.4) 

IF  (N.  LE.  0)  G0T074 
MAXN= ( MAXO ( N, MAXN ) ) 

G ( N, M  >  =GNM 
GT ( N, M ) =GTNM 
GTT ( N, M ) =GTTNM 
IF  (M.  EQ.  1 )  G0T05 
G ( M- 1 , N ) =HNM 
GT ( M— 1 ,  N ) =HTNM 
GTT (M-l,  N ) =HTTNM 
GO  TO  5 
CONTINUE 


c*********WrITE  COEFFICIENTS  FOR  DEBUGGING  PURPOSES 


D7 
D8 
D 

D239 

D 

D 

D 

D 

D 

D9 

D 

DIO 
Dll 
D 1 2 


C*******K 
C*******K 
14 


FORMAT  ?6H0  N  M, 6X1HG,  10X1HH,  9X2HGT, 9X2HHT, 8X3HGTT8X3HHTT// > 
WRITE (5, 239) MAXN 
FORMAT ( 110) 

DO  12  N=2, MAXN 
DO  12  M=l, N 
MI=M— 1 

IF(M.  EQ.  1  )G0  TO  10 
WRITE  (5,9)  N,  M,  G(N,  M) ,  G(MI, 

FORMAT  (213, 6F11.4) 

GO  TO  12 

WRITE  (5,11)  N,  M,  G(N,  M),  GT(N,  M),  GTT ( N,  M ) 

FORMAT  (213,  FI  1.  4,  11X,  Fll.  4,  11XF11.  4) 

CONTINUE 
NE.  0 
EQ.  0 


N  ) ,  GT  (  N,  M ) ,  GT  (MI .  N ) ,  GTT  ( N,  M),  GTT  (MI,  N) 


15 


IF  COEFFICIENTS  ARE  GAUSS  QUASI NORMAL I  ZED 
IF  COEFFICIENTS  ARE  SCHMIDT  QUAS INORMALI ZED 
IF  (K.  NE.  0)  GOTO  17 
SHMIT  ( 1  <_1  >=-!• 

SHM I T  (  n7  1 ')  =SHM  I T  ( N- 1 ,  1  )  *FLOAT  ( 2*N-3  )  / FLOAT  ( N- 1  ) 

SHMIT  (  1 ,  N  )  =0. 

J  J=2 

SHMIT (N, M) =SHMIT ( N, M-l )*SQRT (FLOAT ( (N-M+l ) *JJ ) /FLOAT ( N+M-2 ) ) 
SHMIT (M-l, N) =SHMIT(N,  M) 

JJ=1 

DO  16  N=2, MAXN 
DO  16  M= 1 , N 

G(N, M)=G(N, M)*SHMIT(N, M) 

GT(N, M ) =GT ( N, M)*SHMIT(N,  M> 

GTT i N, M ) =GTT ( N,  M)*SHMIT (N,  M) 

IF  (M.  EQ.  1  )  GOTO  1 6 


TERMS  ARE  DEFINED  AS 


G(M-1, N  >  =G  ( M— 1 1  N ) *SHMIT ( M-l ,  N) 

GT(M-1,  N ) =GT ( M— 1 ,  N ) *SHMIT ( M— 1 ,  N ) 

GTT ( M- 1 ,  N ) =GTT ( M- 1 ,  N  > *SHM I T ( M- 1 ,  N ) 

16  CONTINUE 

17  T=TM— TZERO 
DO  IS  N=1,MAXN 
DO  IS  M=l, N 

TG  (  N,  M  )  =G  (  N,  M  )  +T*  <  GT  (  N,  M  )  +GTT  (  N,  M  )  *T  ) 

CAREFUL  THAT  SECULAR  ACCELERATION 
C **********  A*T/2. 

IF  (M.  EG.  1  )  G0T018 

TG ( M— 1 , N ) =G ( M— 1 , N)+T*(GT(M-1,  N ) +GTT ( M-l ,  N)*T) 

IS  CONTINUE 

TLAST=TM 

19  DLATR=DLAT/57.  2957795 
S I NLA=S I N ( DLATR ) 

RL0NG=DLQNG/57.  2957795 
CPH=COS ( RLQNG ) 

SPH=SIN ( RLONG ) 

IF  (J.  EG.  0)  G0T020 

C ****** !F  J. NE.  0  OUTPUT  IS  IN  GEOCENTRIC  COORDINATES 
R=ALT+6371. 0 
CT=SINLA 
GO  TO  21 

C*****CALCULATE  GEODETIC  CONSTANTS 

20  S I NLA2=S I NLA**2 
C0SLA2=1.  -SINLA2 
DEN2=A2-A2B2*SINLA2 
DEN=SQRT ( DEN2 ) 

FAC= ( ( ( ALT*DEN ) +A2 ) / ( ( ALT*DEN ) +B2 ) ) **2 
CT=SINLA/SQRT ( FAC*C0SLA2+S INLA2 ) 

R=SQRT ( ALT* ( ALT +2.  *DEN )  +  < A4-A4B4*S INLA2 ) /DEN2 ) 

21  ST=SQRT( 1.  — CT**2) 

C******EVALUATE  SPHERICAL  HARMONIC  TO  ORDER  AND  DEGREE  NMAX  (NOTE: 
C******NMAX  IS  N+l  OF  MATH  EXPRESSION  I.  E.  NMAX  14=DEGREE  13  TO 
C*******SATISFY  THE  COMPUTER  DO  LOOPS 
NMAX=MINO(NMX,  MAXN ) 

C*******EVALUATE  SPHERICAL  HARMONICS 
CALL  FIELD 
Y=BP 
F=B 

IF  (J)  22,23.22 
X=— BT 
Z=-BR 
RETURN 

TRANSFORMS  FIELD  TO  GEODETIC  DIRECTIONS 
S I ND=S I NLA*ST-SQRT  <  C  QSLA2 ) *CT 
COSD=SQRT ( 1 .  0— SIND**2> 

X=-BT*COSD-BR*SIND 

Z=BT*SIND-BR*COSD 

D  WR  ITE  (  5,  1212 ) DLAT ,  DLONG,  TM,  X,  Y,  Z,  F 

1212  FORMAT ( 7F8.  1 ) 

RETURN 

END 

SUBROUTINE  FIELD 

COMMON  NDPY ( 3 ) , LPYR ( 4 ) ,  JBSYR 
COMMON/COEFFS/G( 18, 18) , GX ( 18, 18 ) , GT ( 18, 18 ) 

COMMON/FLDCOM/ST , CT, SPH,  CPH, R, NMAX, BT, BP, BR, B 
DIMENSION  P ( 18, 13) , DP ( 18, 18) , CONST (18, 18) , SP ( 18) 

1FM( 18) 

IF  <P(  1,  1  >.  EG.  1.  0)  GO  TO  3 

1  P(  1,  1  )=1. 

DP  (  1,  1  )=0. 

SP  ( 1  )=0. 

CP ( 1 )=1. 

DO  2  N=2, 18 
FN(N)=N 
DO  2  M= 1 , N 
FM ( M ) =M— 1 

2  CONST (N, M)=FLOAT( ( N-2 ) **2- ( M- 1 )**2> /FLOAT ( ( 2*N-3 ) * ( 2*N-5 ) ) 

3  SP ( 2 ) =SPH 
CP ( 2 ) =CPH 
DO  4  M=2, NMAX 

SP (M)=SP (2)*CP (M-l )+CP (2)*SP (M-l ) 

4  CP (M)=CP (2)*CP (M-l )-SP (2)*SP (M-l ) 

A0R=637 1 .  2/R 
AR=A0R*-92 
BT =0. 

BP=0. 

BR=G. 


22 


C 

23 


GTT ( 18, 18), TLAST 


CP ( 18) , FN( 18) 


5 

6 
7 


9 

10 

a 


DO  8  N=2,  NMAX 

AR=AOR*AR 

DO  8  M=1 , N 

IF  (N-M)  6/  5,  6 

P  (N,  N  )  =ST-*P  ( N— 1  /  N— 1  ) 

DP(N,  N ) =ST*DP ( N— 1 1  N-l )+CT*P(N-l, N-l ) 

GO  TO  *7 

P(!M,  M)=CT*P(N-1/  M  )  -CONST  (N,  M)*P(N-2/  M) 
DP  (N,  M ) =CT*DP  (N-l/  M)-ST*P  (N-l/  M) -CONST 
PAR=P(N, M)*AR 
IF  (M.  EQ.  1  )  GO  TO  9 


T£MP=G <  N, M)*CP (M)+G(M— 1/  N)*SP (M> 

BP=BP—  <  G ( N/  M)*SP (M)-G(M-l/  N)*CP(M>  )*FM 
GO  TO  10 

TEMP=G(N,  M)*CP(M> 

BP=BP-(G(N,  M)*SP (M) >*FM(M)*PAR 
BT=BT+TEMP*DP(N, M)*AR 
BR=BR-TEMP-»-FN  ( N  )  *PAR 
BP=BP/ST 

B=SQRT ( BT*BT+BP*BP+BR*BR ) 

RETURN 

END 


(N,  M>*DP(N-2/  M> 

<M)*PAR 
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GSFC 
1000 
5000 
.  5900 
.  9800 
5100 
.  7700 
.  2900 
.  6200 
6740 
.  2870 

0310 

9100 
.  0990 
.  2710 
.  3320 
.  8670 
.  9480 
.  2541 
.  0020 
.  9389 
0175 
.  7832 
.  4079 
.  5830 
.  5753 
.  3706 
.  8067 
.  4081 
.  1911 
.  4095 
.  9282 
.  3203 
.  1347 
9015 
4038 
0202 
0248 
8150 
1829 
5881 
4636 
7581 
8653 
6659 
6246 
1694 
3429 
8288 
9175 
2342 
8658 

3898 
4242 
3130 
2299 
6213 
4385 
2793 
3098 
9974 
3318 
2109 
9023 
7863 
1481 
4064 
2052 
7248 
2204 
1385 
0226 
4953 
6846 
1828 

3899 
4696 
4070 
3164 
5724 
2595 


-1. 

4. 

-3. 


PREMA01  (14,14) 
0.  0 

5834.  5400 
0.  0 

-1981. 6500 
212. 7270 
0.  0 

-438.  3460 
233.  0490 
-174.  5270 
0.  0 

136.  1850 
-278. 7620 
11.  1128 
-241.  4930 
0.  0 

12. 6742 
114.  5650 
-101.  3940 
-106.  6370 
51.  8397 
0.  0 

-10. 2904 
106.  0420 
62.  2257 
-37.  8430 
-19.  0708 
-2.  0597 
0.  0 

-55.  6066 
-26.  8871 
-10.  7983 
3.  1835 
31.  7398 
-2.  5129 
-53.  3257 
0.  0 

7.  3400 
-13.  5012 
7.  3109 
-13. 6024 
14.  5928 
14.  7982 
-6.  2293 

10.  5759 
0.  0 

-22.  3338 
13.  8975 
1. 8808 
-0. 1417 
-9.  4786 
-6.  0963 
22.  7767 

6.  4458 
-34.  5974 

0.  0 

5.  0375 
3.  1052 

1.  5856 

7.  6086 
-4.  1284 

2.  5530 
5.  7351 

-11. 6722 
7.  9897 
9. 4772 
0.  0 

1.  6326 

2.  8587 

2.  4821 
-5.  8210 

3.  9979 
7.  2421 
3.  0473 

-9. 0639 

11.  5000 
0.  6636 

15.  7471 
0.  0 

-0.  5249 
0.  0726 


EPQCH=1960 
23.  7520 
11.  0912 
■23.  6237 
1.  0427 
1.  6995 
-2.  7163 
■10.  9212 
-2.  0714 
-2. 2603 
-1.  1912 
-1.  4004 
-3.  6766 
-1.  3809 
-4.  5806 
1.  0018 
0.  0446 
1 .  3429 
-2.  0854 
-0.  7501 
0.  5085 
0.  3672 
0.  4776 

1 .  7966 

2.  2707 
0.  8038 

1.  5199 
-1. 6451 
-0. 3000 
-0. 0144 

0.  0046 
0.  6028 
1 . 4086 
0.  2407 
-1. 1074 

2.  5005 
0.  1968 
0.  2916 
0.  5765 
0.  8121 

-O. 2094 
-1.  1174 
1 .  2860 
-O.  0513 
-3.  1399 
-0.  0205 
0.  0038 
-O.  4359 
-0.  2164 
-0.  5291 
-0.  9419 
0.  5379 
-0. 3582 
0.  6342 
0.  5547 
O.  0041 
O.  2202 
-0.  1674 
-0. 6292 
-0.  1139 
0.  9014 
-1.  1565 
0.  4386 
0.  7084 
-O.  4505 
2.  3741 
-0. 1746 
0.  0193 
0.  0268 
0.  0571 
0.  5920 
0.  1297 
0.  0235 
0.  3127 
-O. 3453 
0.  4515 
-1.  3535 
-1.  1228 
0.  0168 
-0. 5249 
0.  2161 


.0  A=6341 .  2  F= 

0.  0  0.  0 

-9.  3881  0.  0 

0.  0  0.  0 

-6.  0416  0.  0 

■18.  3544  0.  0 

0.  0  0.  0 

7.  0997  0.  0 

1.  5963  0.  0 

-4.  3583  0.  0 

0.  0  0.  0 

3.  3772  0.  0 

1.0738  0.0 

1.  1670  0.  0 

-2.  8827  0.  0 

0.  0  0.  0 

0.  9806  0.  0 

2.  4208  0.  0 

-3.  0060  0.  0 

2.  0907  0.  0 

1.6224  0.0 

0.  0  0.  0 

-0.  2323  0.  0 

-0.  3307  0.  0 

1.  1895  0.  0 

-O.  4007  0.  0 

0.  3546  0.  0 

3.  5104  0.  0 

0.  O  0.  0 

-1.  3414  0.  0 

-0.  0325  0.  0 

0.  1696  0.  0 

0.  4245  0.  0 

-0.  2976  0.  0 

-1.  6690  0.  0 

0.  4367  0.  0 

0.  0  0.  0 

0.  1595  0.  0 

-0.  1146  0.  0 

-0.  2587  0.  0 

-0.  3901  0.  0 

-0.  4462  0.  0 

-0.  3884  0.  0 

0.  0010  0.  0 

-1.  4777  0.  0 

0.  0  0.  0 

-0.  1244  0.  0 

0.  1891  0.  0 

0.  4797  0.  0 

-0.  1326  0.  0 

0.  1777  0.  0 

1.  3047  0.  0 

-0  2832  0.  0 

-0.  3878  0.  0 

3.  7125  0.  0 

0  0  0.  0 

-0.  3744  0.  0 

-0.  2772  0.  0 

-0.  0597  0.  0 

-0.  1997  0.  0 

-0.  2655  0.  0 

-0.  0234  0.  0 

-0.  6391  0.  0 

1.2878  0.0 

-1.  1128  0.  0 

-0.  5027  0.  0 

0  0  0.  0 

0.  1070  O.  0 

0.  0046  0.  0 

-0.  4512  0.  0 

0.  2162  0.  0 

-0.  1661  0.  0 

-0.  7093  0.  0 

-0.  8235  0.  0 

0.  5943  0.  0 

-1.  1986  0.  0 

0.  2108  0.  0 

-1.  2665  0.  0 

0  0  0.0 

0.  0992  0.  0 

0.  0588  0.  0 
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